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1 Introduction 

The purpose of this article is to explore some simple relations between loop 
measures, spanning trees, determinants, and Gaussian Markov fields. These 
relations are related to Dynkin's isomorphism (cf [1], [11], [7]) . Their poten- 
tial interest could be suggested by noting that loop measures were defined 
in [5] for planar Brownian motion and are related to SLE processes (see also 
[17]). It is also the case for the free field as shown in [13]. We present the 
results in the elementary framework of symmetric Markov chains on a finite 
space, and then indicate how they can be extended to more general Markov 
processes such as the two dimensional Brownian motion. 

2 Symmetric Markov processes on finite spaces 

Notations: Functions on finite (or countable) spaces are often denoted as 
vectors and measures as covectors in coordinates with respect to the canonical 
bases associated with points (the dual base being given by Dirac measures 
6 X ). 

The multiplication operators defined by a function, / acting on functions 
or on measures are in general simply denoted by /, but sometimes multipli- 
cation operators by a function / or a measure Awill be denoted Mf or M\. 
The function obtained as the density of a measure /i with respect to some 
other measure v is simply denoted -. 
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2.1 Energy and Markovian semigroups 

Let us first consider for simplicity the case of a symmetric irreducible Markov 
chain with exponential holding times on a finite space X, with generator 
Ly = q x (Py — 5y), X x , x £ X being a positive measure and P a A-symmetric 
stochastic transition matrix: \ x P y — \P X with P x = for all x in X. 

We denote P t the semigroup exp(Lt) = an d by m x the measure 

. L and P t are m-symmetric. 

Recall that for any complex function z x ,x G X, the "energy" 
e(z) = {-Lz,z) m = -{Lz) x z x m x 

is nonnegative as it can be written 

e( Z ) = ±J2 c ^ zX - zV ^ x - *o + E K * z * zX = E x * z * zX - E 

x,y x a; x, y 

with = = \ X P X and ^ = A x (l - Y. y p y ), i- e - K = «* + Y. y C x, y = 
e(l{*})- 

We say (x, y) is a link iff C x<y > 0. An important exemple is the case of 
a graph: Conductances are equal to zero or one and the conductance matrix 
is the incidence matrix of the graph. 

The (complex) Dirichlet space H is the space of complex functions 
equipped with the energy scalar product defined by polarisation of e. Note 
that the non negative symmetric "conductance matrix" C and the non nega- 
tive equilibrium or "killing" measure k are the free parameters of the model, 
(so is q but we will see it is irrelevant for our purpose and we will mostly 
take it equal to 1). The lowest eigenvector of — L is nonnegative by the well 
known argument which shows that the modulus contraction z — > \z\ lowers 
the energy. We will assume (although it is not always necessary) the corre- 
sponding eigenvalue is positive which means there is a "mass gap": For some 
positive e, the energy e(z) dominates e (z,z) rn for all z. 

We denote by V the associated potential operator (— L)^ 1 = J °° P t dt. 
They can be expressed in terms of the spectral resolution of L. 

We denote by G the Green function defined on X 2 as G x ' y — — 

J m y 

- Py l Y y i.e. G — (M\ — C)-\ It verifies e(f, Gfj,) = (/,//) for all 
function / and measure /i. In particular Gk = 1. 
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Different Markov chains associated to the same energy are equivalent un- 
der time change. If g is a positive function on X, in the new time scale 
f*g£ s ds, we obtain a Markov chain with gm-symmetric generator ^L. Ob- 
jects invariant under time change are called intrinsic. The energy e, P and 
the Green function G are obviously intrinsic but L, V and P t are not. We 
will be interested only in intrinsic objects. In this elementary framework, it 
is possible to define a natural canonical time scale by taking q — 1, but it 
will not be true on continuous spaces. 

2.2 Recurrent chain 

Assume for simplicity that q = 1. It will be convenient to add a cemetery 
point A to X, and extend C, A and G to X A = {XU A}by setting C x> a = k x 
> A a = E, G x «*■ and G X > A = 0. Note that A(X A ) = £xxx C x , v + 2 £ x k x ) 

One can consider the recurrent " resurrected" Markov chain defined by the 
extensions the conductances to X A . An energy e R is defined by the formula 

e R {z)= l -Y, C ^ zX - zV )^-^) 

We denote by P R the transition kernel on X A defined by 

e R (z) = (z-P R z,z) x 

or equivalently by 

tpRlX _ C X ,y _ C X) y 

Note that P K 1 = 1 so that A is now an invariant measure. Let A -1 be the space 
of functions on X A of zero A measure and by V R the inverse of the restriction 
of I — P R to A -1 . It vanishes on constants and has a mass gap on A -1 . Setting 
for any signed measure v of total charge zero G R v = V R j. we have for any 
function /, (u, f) = e R (G R u, f) and in particular f x -p = e R (G R (5 x -5 y ), f). 

Note that for \x E X 1 - and carried by X, for all x G X, /i x = e R (G R fi, l x ) = 
X X ((I — P)G R a){x) — K x G R fi(A). Hence, applying G , it follows that on 
X, G R fi = G r u,(A)Gk + G/i = G r /j,(A) + Gfi. Moreover, as G R fi is in A x , 
G R fx(A)X(X A ) + Y^xex ^x(Gfi) x = 0. 

Therefore, G R /i(A) = =MQ and G R fi = ^g^f + Gfi 
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2.3 Transfer matrix 

We can define a scalar product on the space A of antisymmetric functions 
on X A x X A as follows 

(uj,rj) = J2xyCx,yU x ' v ri x ' y . Denoting as in [9] df u ' v — f u — f v , we note 
that (df,dg) = e R (f,g) In particular 

(df,dG R (5 x -5 y ))=df x >y 

As the antisymmetric functions df span the space of antisymmetric functions, 
it follows that the scalar product is positive definite. 

The symmetric transfer matrix K, indexed by pairs of oriented links, is 
defined to be 

K (x,y),(u,v) = G R^ x _ 5y) u _ qR^ _ Sy y =< dG R( S x _ gy^ dQ R^u _ ^ > 

for x, y,u,v G X A , with x ^ y, u ^ v . 

We see that for x and y in X, G R (5 X - 5 y ) u - G R (5 X - S y ) v = G(5 X - 
5 y ) u - G(5 X - 5 y f. 

We can see also that G R (5 X — 5a) = GS X — ~j^x^f • ^° ^ ne same identity 
holds in X A . 

Therefore, as G X,A = 0, in all cases, 

J^(x,y),(u,v) _ QX,u _|_ Qy,v _ qx,v _ Qy,u 

For every oriented link £ = (x,y) in X A ,set = dG R (5 x — 5 y ) = 
dG(5 x -5y). 

We have (K^, K v } = K^ ,v . K will be viewed as a linear operator on A, 
self adjoint with respect to (•, •) . (It can also be viewed as symmetric with 
respect to the euclidean scalar product if we wish to use it Then it appears 
as the inverse of the operator defined by (•,•)). 

3 Loop measures 
3.1 Definitions 

For any integer k, let us define a based loop with p points in X as a couple 
(f , r) = ((f m , 1 < m < p), {r m , 1 < m < p + 1), ) in X p x and set 

£i = £ p +i- p will be denoted 



4 



Based loops have a natural time parametrisation and a time period 

T(0 = Effi +1 ^- If we denote YZi T * ^ T m . £(t) = U-i on [T m . u T m ) 
(with by convention T = and £o — £p)- 

A cr-finite measure /x is defined on based loops by 



= E / 7^ ,Xdt 



where F^ ,x denotes the (non normalized) "law" of a path from x to x of 
duration t : If Yli=i ^ = ^ 



a; 



pr(^i) = x 1 ,...,t(t H ) = x h ) = [p^uPt 2 -^i\..\Pt-t h . 

Note also that 

K X {P — k,^2 — x 2 , £ fc = x k , T x G tftx, T fc G dt k ) 

= [P]:JP]^^..[P]^n {0<^1< ... tfc<t} g^-^ ^ ^..g Xfc e-^ fc ^- ^ -^e-^^ t ^rf^ 1 ...^ 

A loop is defined as an equivalence class of based loops for the R-shift that 
acts naturally. u. is shift invariant, It induces a measure \i on loops. 

Note also that the measure du = l ,J' qil — du which is not shift invariant 

Jo 1tW ds 

also induces [i on loops. 
It writes 



A*o(p(0 = k,^ = Xi, ...,£ k = x k , Ti G dt ± , T fc G dt k , T G dt) 

0<ti<...<t fc . 



for /c > 2 and 



A*o{p(0 = 1,6 = a:,Ti e <&i} = — : — ^i 

It is clear, in that form, that a time change transforms the /Vs of Markov 
chains associated with the same energy one into each other, and therefore 
the same holds for /j,: this is analogous to conformal invariance. Hence the 
restriction /// of \i to the a-field of sets of loops invariant by time change (i.e. 
intrinsic sets) is intrinsic. It depends only on e. As we are interested in the 
restriction \ij of ji to intrinsic sets, from now on we will denote simply /ij by 
fi 
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Intrinsic sets are defined by the discrete loop £j (in circular order, up to 
translation) and the associated intrinsic times — = r*. Conditionally to the 
discrete loop, these are independent exponential variables with parameters 

M = E^- T '^ + E E II C^e-^dr- (1) 

Sets of discrete loop are the most important intrinsic sets, though we 
will see that to establish a connection with Gaussian fields it is important to 
consider occupation times. The simplest intrinsic variables are 

N x , y = #{i : & = x,£ i+1 = y} 

and 

N X = ^N X ,y 

y 

Note that N x = > 1 : & = x} except for trivial one point loops. 

A bridge measure /i x ' y can be defined on paths 7 from x to y: n x,y (d^) = 

A-r-p-(d 7) ^with 



^t' y (^(ti) =x 1 ,...,'y(t h ) =x h ) = P tl (x,x 1 )P t2 - tl (x 1 ,x 2 )...Pt-t h (xh,y) 



v x 

Note that the mass of u x,y is — = G x,y . We also have, with similar notations 

m y 

as the one defined for loops 



H x,y (p(l) = K 72 = x 2 , 7fc_i = Ti G dti, T fc _! G cft fe _i, T G dt) 
r r n 

^ X ,X2^ X2,X3 - ■ - ^ Xk ^ 1 ,y . 



X X A X2 ...Ay 



-l{o< tl <...< tfc <t}e- fe ' 1 e-^ (i2 - il) ...e- <? " (t - tfc) g :E rft 1 ...^ fc _ 1 rft fc g^t 



so that the restriction of fj, x,y to intrinsic sets of paths is intrinsic. 
Finally, we denote the family of probability laws on paths defined by 

Pt- 

^ x {l{t\) =x 1 ,...,'y(t h ) =x h ) = P tl (x,x 1 )Pt a -t 1 (x 1 ,x 2 )...Pt h -t h _ 1 (x h - 1 ,x h ) 
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...e-^ tk - t ^qxdt 1 ...q Xk dt k 



3.2 First properties 

If D is a subset of X, the restriction of \x to loops contained in D, denoted 
fi D is clearly the loop measure induced by the Markov chain killed at the exit 
of D. This can be called the restriction property. 

Let us recall that this killed Markov chain is defined by the restriction of 
A to D and the restriction P D of P to D 2 (or equivalently by the restriction 
eo of the Dirichlet norm e to functions vanishing outside D) and (for the 
time scale), by the restriction of q to D. 

From now on in this section, we will take q x — 1 for all x. Then /x takes 
a simpler form: 



for fc > 1 and /U {p(£) = 1, £i — x±,Ti E dt\\ = ^-dt\ 
It follows that for k > 0, 



'.T 



-AT, 



as/ t—Q-tdt — I anc [ conditionally to = = xi,...,£fc = T is a 
gamma variable of density (fz^y e - * on M + and (^-1 < i < k) an independent 
ordered k— sample of the uniform distribution on (0, 1). 



In particular, we obtain that, for k > 2 



/*(p — k) — fi (p = k) = l -Tr{P k ) 



and therefore, as Tr{P) = 0, 



ji(p > 0) = - log(det(7 - P)) = - log( TT V) 
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as denoting M x the diagonal matrix with entries \ x , det(I-P) = dc ^^ f ^^ . 
Moreover 



jp(l)f,(dl)=Tr((I-P)- 1 P) 



Similarly, for any x ^ y in X and s G [0, 1], setting Pu} = P™ if (u, v) ^ 
(x,y) and P x ^ = sP x , we have: 

M^l {p>0} ) = -log(det(J - ?W)) 
Differentiating in s = 1, it comes that 

/<(.V„„) = [(/ " P)- l ]lP X y = G^C X ,y 

and »(N X ) = J2 y K N x, y ) = KG X ' X - 1 (as (M A - C)G = Id). 

4 Poisson process of loops and occupation 
field 

4.1 Occupation field 

To each loop I we associate an occupation field {l x ,x G X} defined by 

f T{l) ^ q x n ^ 

F = 1 {as)=x}——ds = 2^ %_!=*} — ~ = 2^ l {Hi-i=x} T i 

Jo m i(s) m x , =1 



for any representative (£, t) of/. It is independent of the time scale (i.e." intrinsic"). 
For a path 7, 7 is defined in the same way. 
From now on we will take q — 1. 
Note that 

M(l - e-nW = /°°(e-^^ - e"*)^ = M^-) (2) 
In particular, /i(/ rr l{ p=1 }) = j-. 

From formula 1, we get easily that for any function $ of the discrete loop 
and k > 1, 

M(n fc l {p >i}$) = A*((JV* + k - 1)...(N X + 1)N X $) 



In particular, nff x ) = j^\fi(N x ) + 1] = G x ' x . 

Note that functions of I are not the only intrinsic functions. Other in- 
trinsic variables of interest are, for k > 2 

yc\,...,x k _ _ £\ =Q J 0<tl< <tk<T l{S(ti)=x 1+j ,....S(t k - j )=x k ,...€(t k )=x j } II A^7^» 

= I EjU Ei<i 1 <..<i fc < P (i) nti l{? i( _i=WT* and one can check that fitf* 1 --**) = 
G Xl ' X2 G X2,Xs ...G Xk ' Xl . Note that in general l x i>-> x k cannot be expressed in 
terms of / for k > 3. 

For xi = x 2 = ... = Xk, we obtain self intersection local times l x ' k = 

^2l<h<..<i k <p(l) Yll=l ^{&[-i=x} r i t 

For any function $ of the discrete loop, fj,(l x ' 2 &) = //( ^ £ ^— ^ <fr) since 
p, 2 = I((?) 2 - E g) l te _ 1=a} (r;) 2 ) and/,($ES Ife-^*) 2 )) = 2^($JV X ) 

More generally one proves in a similar way that n{l x,k <&) = ^{ Nx<yNx ^'^ Nx fc+1 ^ «l>) 

From the Feynman-Kac formula, it comes easily that, denoting Mi the 
diagonal matrix with coefficients ^ 

K,x(e-( T ' X ) -1) = exp(*(P-J-M $ )) X)X -exp(*(P-J)) X)X . Integrating in 
t after expanding, we get from the definition of \x (first for x small enough): 

/ OO -i 
(e-(^) - = g -[Tr((P - Mj fc ) - Tr((P) fc )] 

Hence 

/(e-<-> - 1 )<W0 = .o g [de t (-L(-L + M^)-)] = - log det(/ + KM,) 

which now holds for all non negative x- Set V x = (-L + M x ) _1 and G x = 

V x Mi. It is an intrinsic symmetric nonnegative function on X x X. Gq is 
the dreen function G, and G x can be viewed as the Green function of the 
energy form e x = e + \\\\ 2 L 2/ X y Note that e x has the same conductances C 
as e, but x is added to the killing measure. We have also the "resolvent" 
equation V -V x = VMxV x = V x MxV. Then, G-G x = GM X G X = G X M X G. 
Also: 

det(J + GM X ) _1 = det(J - G X M X ) = ^2*) (3) 
Finally we have the 
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Proposition 1 i>(e"^ x )-l) = - log(det (J+GM X )) = log(det( J-G X M X )) = 
log(det(G x G- 1 )) 

Note that in this calculation, the trace and the determinant are applied to 
matrices indexed by X. Note also that det(J + GM X ) = det(J + M ^GM ^) 
and det(J — G X M X ) = det(I — M ^ X G X M so we can deal with symmetric 
matrices.. 

In view of generalizing them to continuous spaces in an intrinsic form (i.e. 
in a form invariant under time change), , G and G x will be interpreted as 
symmetric elements of H <g) H, or as linear operators from H 7 into H. G is a 
canonical bijection. ^r^j can be viewed as the determinant of the operator 
G X G _1 acting on H. 

4.2 Poisson process of loops 

Still following the idea of [5], define, for all positive a, the Poisson process 
of loops C a with intensity aji. We denote by P or F Ca its distribution. Note 
that by the restriction property, C% = V e ^ ^} i s a Poisson process 
of loops with intensity fj, D , and that is independent of C a \C^. 

We denote by £^ the set of non trivial discrete loops in C a . Then, 

na = {h,k,...i k }) = e -^a k ^^ = [ggriLXr n x K N * a) 

with # = J2 leCa Nx (l) and iVg = ^c a N ^)- 

Remark 2 It follows that the probability of a discrete loop configuration de- 
pends only on the variables N XjV + N y ^ x , i.e. the total number of traversals 
of non oriented links. In particular, it does not depend on the orientation of 
the loops It should be noted that under loop or path measures, the conditional 
distributions of discrete loops or paths given the values of all N x ^ y + N y>x 's 
is uniform. The N x ^ y + N y>x (N x>y ) configuration can be called the associated 
random (oriented) graph. Note however that any configuration of N x ^ y + N y ^ x 
does not correspond to a loop configuration. 

We can associate to C a the a-finite measure 

Then, for any non-negative measure x on X 

E( e -<£'*)) = exp(« / (e-( T > x ) - l)dfi(l)) 
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and 

E ( e -(^,x)) = [det(-L(-L + M^a)- 1 )]" = det(J + VMx)~ a 
Finally we have the 

Proposition 3 E(e - ^'*>) = det(/+GM x )- a = det(/-G x M x ) a = det(G x G'- 1 ) a 

Many calculations follow from proposition 1. 

It follows that K(C a ) = aG X)X and we recover that fi(l x ) = G XyX . 

On loops and paths, we define the restricted intrinsic a-field Xr as gener- 
ated the variables N x>y with y. possibly equal to A in the case of paths, with 
N xA = or 1. from (2), 

k 

The distribution of {N x a \x G X} follows easily, in terms of generating 
functions: 



i=i V s<6 j 

Note also that 

FrrrV|T , (iVi a) + A:)(iVi a) + A:-l)...(iVi") + l) 

mCa } |Jr) = m 



andjf se_lf intersection local times are defined as 

= ELi E fcl +...+fc m =fc E/ 1 ^ 2 ...^i m e£+ nr=i J ' we § et easil y that 



E(C a X ' k \l R ) = 1(JVW - fc + 1)...(JV<°> - 1)7V>) 



Note also that since G X M X is a contraction, from determinant expansions 
given in [15] and [16], we have 

) = ^2xh--Xi k Per a (G ihim ,l <l,m<k) 
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Here the a-permanent Per a is denned as Xlo-G5 fc am ^^u,V(i)--^»k>Mfc) 
with m(o~) denoting the number of cycles in a. 

Let [H F \ X be the hitting distribution of F by the Markov chain starting 
at F. Set D = F C and denote e D , V D = [(/ - P^dxd}' 1 and G D = [(M A - 
C)\d x dY 1 the Dirichlet norm, the potential and the Green function of the 
process killed at the hitting of F. Recall that V = V D + H F V and G = 
G D + H F G. 

Taking x — a ^F with F finite, and letting a increase to infinity, we get 
lim a | 00 (G x M x ) = H F which is / on F. Therefore by proposition 1, one 

checks that F{C a {F) = 0) = det(J - H F ) = and n(T(F) > 0) = oo. But 
this is clearly due to trivial loops as it can be seen directly from the definition 
of fj, that in this simple framework they cover the whole space X. 

Note however that n(t(F) > 0,p > 0) = fi(p > 0) - n(t(F) = 0,p>0) 

= fi(p > o) - fip > o) = - hg( = M njSM ) 

It follows that the probability no non trivial loop (i.e. a loop which is not 
reduced to a point) in C a intersects F equals ( ^ ^^detfG) )" 

Recall that for any (n + p, n+p) invertible matrix A, det(A _1 ) det(Ajjl < 
i,j <n) = det(A~ 1 ) det(Aei, ...Ae n , e n+1 , ...e n+p ) 

= det(ei, ...e n , A _1 e n+ i, ...A _1 e n+P ) = det((^4 _1 ) fc) j, n < k,l <n + p). 

In particular, det(G D ) = dc ^Q\^} F y so we have the 

Corollary 4 The probability that no non trivial loop in C a intersects F 
equals (H xeF \ x det FxF (Gy a 

In particular, it follows that the probability a non trivial loop in C a visits 
x equals 1 - ( Xx Q X , x ) a 

Also, if Fi and F 2 are disjoint, /x(f] l(FA > 0) = fi(p > 0) + / (Fj) = 
0,p > 0) - n(T(Fi) = O p > 0) - //(f(F 2 ) = 0,p > 0) 

log( ^(g^i ) dct(G°2^ ) an d this formula is easily generalized to n disjoint 

- det(G) n^, det(G^ nD 0... 

Mil W > °) = l0g( ndet(^)n, <J<fc det(^.n Dfc) ^ 

The positivity yields an interesting determinant product inequality. 

It follows in particular that the probability a non trivial loop in C a visits 
two distinct points x and y equals 1 — ( G ' qI'^q^I ) a and q^ x 'q v , v if a = 1. 



sets. 
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Note finally that if x has support in D, by the restriction property 



Ml { ,»=o } ( e ' X> -!)) = - Iog(det(/ + G m x )) = log(det(G^)[G D ]- 1 ) 

Here the determinants are taken on matrices indexed by D. or equiva- 
lently on operators on H D . 

For paths we have P^(e~( i,x )) = exp(t(L - M x )) x>y . 
Hence /^( e -<™>) = JL((j _ p + M^) -1 )^ = [G x ]*.». 
Also E*( e -<™>) = 

In the case of a lattice, one can consider a Poisson process of loops with 
intensity fif 



5 Associated Gaussian field 

By a well known calculation, if X is finite, for any x £ 

det(M A -C) r-.i <rty> ,_i eWn A ^ det(G, 



y ( e -§<«.x> e -|ewn ttex id^ A dz u 



(2tt)w y v ~ u ^ 2 ^ u ™„ det(G) 

and 

det(M A + M x -C) / , ^ l/ ._i<- Y> ._i ef ^ TT < 



(2tt)I x I 



J z x z y (e-^< zz ^> e -^Il ueX ^dz u A cfe" = (G x ) x ' y 



This can be easily reformulated by introducing the complex Gaussian field 
4> defined by the covariance E^(0 a: y ) = 2G x,y (this reformulation cannot be 
dispensed with when X becomes infinite) 

So we have E((e-3<^>*>) = det(J + GMJ" 1 = det^G" 1 ) and 
E(((^0V5«^x>) = (G x ) x » det^G" 1 ) Then the following holds: 

Theorem 5 a) The fields C\ and \<\>§ have the same distribution. 

b) E+((<j) x lj? F(<j>lj))) = / E(F(£ + 7))/j, x ' y (d-f) for any functional F of a 
non negative field. 

This is a version of Dynkin's isomorphism (Cf [1]). It can be extended to 
non symmetric generators (Cf [10]). 
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Note it implies immediately that the process (jxfi is infinitely divisible. See 
[2] and its references for a converse and earlier proofs of this last fact. 

In fact an analogous result can be given when a is any positive half 
integer, by using a real scalar or vector valued Gaussian field. 

Recall that for any / el, the law of / + is absolutely continuous with 
respect to the law of 0, with density exp(< —Lf,<p > m — |e(/)) 

Recall (it was observed by Nelson in the context of the free field) that 
the Gaussian field <fi is Markovian: Given any subset F of X, denote Tip 
the Gaussian space spanned by {(fi y ,y € F}. Then, for x G D = F c , the 
projection of (j) x on Hf is Y^yeF^^y^ ■ 

Moreover, 4> D = 4> — H F 4> is the Gaussian field associated with the process 
killed at the exit of D. 

Note also that if a function h is such that Lh < 0, the loop measure 
defined by the /i 2 m-symmetric generator L h = \LMh is associated with the 
Gaussian field h(j). The killing measure becomes 

Remark finally that the transfer matrix K is the covariance matrix of the 
Gaussian field d(f) x ' y = (j) x — 4> y indexed by oriented links. 



6 Energy variation and currents 

The loop measure \i depends on the energy e which is defined by the free 
parameters C, k. It will sometimes be denoted \i e . We shall denote Z e the 
determinant det(G) = det(M A -C) -1 . Then /x(p > 0) = log(^ e )+X)log(A :r ). 

Other intrinsic variables of interest on the loop space are associated with 
real antisymmetric matrices oj x>y indexed by X A : u XtV = —uj y , x .. Let us 
mention a few elementary results. 

The operator [P u ] x y = P x exp(iu X:V ) is self adjoint in L 2 (A).The associated 
loop variable writes Y7 j= Y^x,y u x,yN x , y {i). We will denote it j\uo. 

This notation will be used even when uo is not antisymmetric. Note it is 
invariant if u x>y is replaced by uj x ^ y + g{x) — g(y) for some g. Set [G fa, ] x ' s/ = 

-r and denote Z e ^ the determinant det(G a; ). By an argument similar 

to the one given above for the occupation field, we have: 

KM SlU ~ X ) = e Mt( puJ ~ - exp(t(P - I)) XjX . Integrating in t 

after expanding, we get from the definition of \i : 

(jfi" - l)dn(l) = J2Tl T <(n k )-Tr((P) k )] 



14 



Hence 

J (e^ w - l)dfx(l) = log[det(-L(7 - P^y 1 } 

and 

KeMYl i I") ~ l ) = log(det(G^G- 1 )) = log(%^) (4) 

The following result is suggested by an analogy with quantum field theory 
(Cf[3]). 

Proposition 6 i)-§£- = T x fi 

™) d\ogC x ,y = ~T x ,yli 

With T Xt y(l) = C X ,y(?+fy) - N x ,y(l) - Ny, X (l) 

Note that the formula i) would a direct consequence of the Dynkin iso- 
morphism if we considered only sets defined by the occupation field. 

Recall that // = £ f x e~ x ^ £+E~ 2 £(^ e z/ P z) e x, FUz/pZ C^ i+1 e~ x ^ 
C x ,y = Cy )X = ^ x Py and \ x = k x + C Xt y 
The formulas follow by elementary calculation. 
Recall that fiff x ) = G x ' x .and n(N X:y ) = G x > y C x , y 
So we have fi(T x>y ) = C x>y (G x > x + G y ' y - 2G x ' y j 

Then, the above proposition allows to compute all moments of T and I 
relative to /i e (Schwinger functions) 

Consider now another energy form e' defining an equivalent norm on H. 
Then we have the following identity: 

d/V = e Eiv,, y iog(g^)-E(A^-A ;c )P 

The above proposition is the infinitesimal form of this formula. Note that 
from the above expression of /i (??), 

/ie((e E^.lo g (if)-E(A-A,)P _ 1)} = , og( ^ ) 

A3 

(the proof goes by evaluating separately the contribution of trivial loops, 
which equals XLlog(ff))- 

Note that if C xy = h x h y C x>y et k' x = for some positive function h 

on E such that Lh < 0, ^ = y^K*)* ■ 
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Note also thatff = E(e-^ e '- e ^) 
Equivalent ly 



(x,y) X ' y x A * x,y r v x ^ 



and therefore 



(5) 



(*,y) Cx ' y x x Ze 

Note also that U {x , y) [^] N ^ = U {x , y} [^} N ^ +N ^ 

N.B.: These ^ determine, when e' varies with ^=r < 1 and — = 1, the 
Laplace transform of the distribution of the traversal numbers of non oriented 
links N X;V + N y;X , hence the loop distribution /i e . 

More generally 

/ie(e -E^io g (^)-E(AL-A, ) r, +i v _ 1} = log( ^v } (6) 



or 



Note also that this last formula applies to the calculation of loop indices if 
we have for exemple a simple random walk on an oriented two dimensional 
lattice. In such cases, u z i can be chosen such that Jj u z > is the winding 
number of the loop around a given point z' of the dual lattice 1 X'. Then 
e z n J2 l G C a J. uj' z is a spin system of interest. 
We then get for exemple that 



H( !u) ^ 0) = -— / log(det(G 2 ^G'- 1 ))rfM 

./z 27T J 



1 The construction of u> can be done as follows: Let P' be the uniform Markov transition 
probability on neighbouring points of the dual lattice and let h be a function such that 
P'h = h except in z'. Then if the link xy in X intersects x'y' in X', with det(x—y, x'—y') > 
0, set u XlV = h(y') - h(x') 
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and hence 



P( I fu'z\) = Q ) = e^Mdct(G 2 — G-i))du 
iec a Jl 



Conditional distributions of the occupation field with respect to values of 
the winding number can also be obtained. 

We can apply the formula 5 to calculations concerning the links visited 
by the loops (similar to those done in section 4 for sites). 

For exemple, R is a set of links, denote e' R [ the energy form defined from 
e by setting all conductances in R to zero and increasing k in such a way 
that A is unchanged.. 

Then He{Y,{ x , y )eR N x,y + N w > 0) = - log( d ^ e ( t ^ 1) ) and therefore, the 
probability no loop in C a visits R equals de ^^ G ^ = (~z^) a - 

7 Self-avoiding paths and spanning trees. 

Recall that link / is a pair of points (f + ,f~) such that C/ = C/+ /- ^ 0. 
Define -/ = (/",/+). 

Let /4f„ be the measure induced by C on discrete self-avoiding paths. between 
x and y: fj,_l y (x : X2, Xn-i, y) = C x>X2 C xl:X3 ...C Xn _ 1:y . 

Another way to defined a measure on discrete self avoiding paths from x 
to y is loop erasure (see for exemple [4]). One checks easily the following: 

Proposition 7 the image of fi x ' y by the loop erasure map 7 — > r y BE is fi x EE 
defined on self avoiding paths by ^ x B %(r]) = fJ^ y (v) d ^St c ) = ^/(v) det(G|{„} X {i,}) 
(Here {r]} denotes the set of points in the path rj) 

Proof: If 77 = (xi = x,x 2 ,...x n = y),and r] m = (x,...x m ), [i x > y (-f BE = 

\V{Vrn-l} C ]Xm = det([(J-P] | {t?m} c x {??m }e ) dct( y{^-l} C ) dct } C ) m 

L V Js m det([(/-P]| {t)m _ l} c x{l)m _ l}c ) - det(V{^} c ) " det(G{W c ) 

all m < n — 1. 

Also:/ e-<^>l {TSB=r)}/ ^(^) = -^L e -^>^/( V ) 



det(G™ ) 

= det(G x ) |WxW e-<^> /i ^(^) = g^gffi e-<^>^| (?7 ) for any 
self-avoiding path rj. 
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Therefore, under fj, x ' y , the conditional distribution of 7 — rj given r y BE = r\ 

is the distribution of L\ — C\ v ^ c i.e. the occupation field of the loops of Ci 
which intersect 77. 

More generally, it can be shown that 

Proposition 8 the conditional distribution of the set of loops of 7 given 
r y BE = rj is the distribution of Ci/C[ v ^ i.e. the loops of C\ which intersect 
V- 

Proof: First an elementary calculation shows that 

»7b BE = v) = g^Sii^Siiii ^nn^j^]^-^^^-^^ m^r^ 

Therefore, by the previous proposition, 

u x,y(Y\ \9^2L\N x , y {C)+N y , x {C) TT \\v/\N u (C)\ n ,BE _ \ _ ZeZ [e']M c 

Moreover, by 5 and the properties of the Poisson processes, 

Kf FT \^]^,y{Ci/c^ }C )+N y , x {c 1 /c\ ll}C ) tt rA^iJVu^i//:^) _ f^W^ 

It follows that the distributions of the N Xjy +N yj xS are identical for the set 
of erased loops and C\/C^ . Moreover, remark 2 allows to conclude, since 
the same conditional equidistribution property holds for the configurations 
of erased loops. 

Similarly one can define the image of F x by BE which is given by 
^be(v) = C XuX2 ...C Xn _ ltXn K Xn det(G\ {v}x{r)} ), for rj = {x u x n ), and get 
the same results. 

Wilson's algorithm (see [9]) iterates this construction, starting with x's 
in arbitrary order. Each step of the algorithm reproduces the first step 
except it stops when it hits the already constructed tree of self avoiding 
paths. It provides a construction of the probability measure P| r on the set 
STx,a of spanning trees of X rooted at the cemetery point A defined by the 
energy e. The weight attached to each oriented link £ = (x, y) of X x X 
is the conductance and the weight attached to the link (x, A) is k x . As 
the determinants simplify, the probability of a tree T is given by the simple 
formula 

Psr(T) = Z e J] C € 
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Proposition 9 The random set of discrete loops C\y constructed in this al- 
gorithm is independent of the random spanning tree, and independent of the 
ordering. It has the same Poisson distribution as the non trivial discrete 
loops of C\. 



It follows easily from proposition 8 . 

Together with the spanning tree these discrete loops define an interesting 
random graph. 

First note that, since we get a probability 



or equivalently 



Z e * ~ j I C Xy y || k x — 1 
?eST X:A (x, y )er x,(x,A)er 



e n >? n ^n^u 

r&ST x ,A(x,y)£T x,(x,A)eT 1 lx ^ X 



so that 

n T m = Zei[K n p v n p a 

xex (x, y )er x,(x,A)eT 
Then, it comes that, for any e', 

p'x p/x TT \ 7 

E | T ( tt — TT — ) = lr gx — 

STy 11 px 11 px > TT \/ 7 , 

(x,»)6T x,(x,A)6T A 1 l*eX A x <=-e' 

and 

»« n I; n J) = | m 

(x,j/)eT ' y x,(x,A)6T 

We also have F e ST ((x,y) E T) = Pf B (^i = y) = V£P£F«(T X = oo) = 

C x ,yG (1 — 

From the results exposed in [8] and [9], or directly from the above, we 
recover Kirchhoff's theorem: 

P| r (±(x, y ) g T) = r,,, y Y,-d - g£) + g**(i - = r,,, /( r;- + 

G y ' y — 2G x ' y ) = C X: yK x ' y ^' ( - x ' y ^ and more generally Pemantle's transfer current 
theorem: 
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k 

nAHi,- ±& e T) = ([[C^det(K^ 1 < i,j < k) 

1 

Note this determinant does not depend on the orientation of the links. 

Proof: We use recurrence on k. Let M denote the smallest subset of 
X A containing the links ±£i, ... ± and denote E — M by D Let V be the 
subspace of A spanned by all K^'^ with x and y in M k . Note that the 
orthogonal of V in A is spanned by dG D (S x ) and that for any rj = (u, v) the 
projection of K" on V L is dG D (6 v -6 u ), and (dG D (S v - 5 u ),dG D (5 v - 5 U )) = 
[QDpu + [QDy,,v _ 2[G D ] U ' V 

Moreover det(K^ 1 < i,j < k) = \\K^ A ... A K^ k ||^ fcA - Therefore, if 
£ fc+1 = n det(K&& 1 < i, j < k + 1) = det(f^ 1 < i, j < k)([G D ] u ' u + 
[G D ] V ' V — 2[G D ] U ' V ) 

But the argument given for k — 1 shows also that P| T (±r] e T| ±£i, ... ± 
£ k e T) = C U ,„([G D ]"' M + [G D ] V ' V - 2[G D } U ' V ) so we can conclude. 

Therefore, given any function g on non oriented links, E| T (e~^-'« eT9 ^- ) ) = 
E e 5 T (n 5 (l + (e~ 9(5) - 1)1 C 6t) = E^((^c(e- 9 -i)^) Afc ) and we have 

E| T ( e -^^(0) = det(/+KM c(e - 9 _ 1) ) = det^-il^— -KM^^_ -) 

This is an exemple of the Fermi point processes discussed in [14]. 
But, by (7) and (5), it comes that 



io g (Ei T ( n 5* n ^)=-Mn[^]^n[#^ +i - i )=M^) 



(x,y)&T x ' y x,(x,A)eT 1 x,y 



The first identity could also be derived from proposition. 9. 

As #■ = ft r, (with the convention — = A~ if « x = 0) we 

obtain 

For any function (7 on non oriented link of Xa, non negative on links of 

X 

E Ca (e-Z*,y9«*'y}) N ™ JJ[ P*e- 9{x > z} ]- N *- r ) = det(J - KM c[1 _ e - g) )- a 

x zex A 
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We can check that this formula allows to recover the identity E^ a ( - A v 

N x , y - N y , x ) = aC x , y (G x > x + G y ' y - 2G x ' y ) It also gives back prpositin 3 for 
g = onX x X. 

If k is positive everywhere, we can adjust g({x, A}) to make J2 zeX A P x e^ 9 ^ x 
1 This means we have to choose re x (l-e- fl « x ' A ») = « x (l-^(l-X) zeX P x e~ 9 ^ x 
= (^-X x +-E z exC xz e-^y) = j: zeX C xz (e- y ^ 
We check also also that by 5 

Eca ( e -i: x , y 9({x,y})N x , y j = = ( de ^-^ 9 ) )-« Finally the restriction 

on k can be removed by taking a limit and we obtain: 

Proposition 10 For any function g on non oriented link of non neg- 
ative on links of X , set Tg(£) = Q(l — e~ 9< ^) if £ is a link of X and 
Tg({x, A}) = Y, z ex C xz (e~^ x ^ - 1) for all x. Then 

det(J + C(I - [e- g ]))- a = E Ca (e-Z«>v 9 U x > v » N *») = det(J - KM(g)y a 

We see that the Poisson measure on loops C a induces a point process iV 
on the space of non oriented links defined by the pair (a, K) which reminds 
the point processes discussed in [14]. Note however a difference of sign in the 
right hand side determinant, which is not a Laplace transform for positive a. 

8 Fock spaces and Wick product 

Recall that the Gaussian space H. spanned by {<p x , x e X} is isomorphic to EI 
by the linear map mapping R,e(<f) x ) on G Xj . which extends into an isomorphism 
between the space of square integrable functionals of the Gaussian fields 
and the symmetric Fock space obtained as the closure of the sum of all 
symmetric tensor powers of H (Bose second quantization). We have seen 
that L? functionals of C\ can be represented in this symmetric Fock space. 

In order to prepare the extension of these isomorphisms to a more interest- 
ing framework (including especially the planar Brownian motion considered 
in [5]) we shall introduce the renormalized (or Wick) powers of 00. 

The Laguerre polynomials L° n (x) are defined by their generating function 

^ = Et n L° n (x) 

Then one defines the polynomial P n {-) = (— l)"n!L° (■) 
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Setting a x = G x , x , it comes that cr n P n (^-) is the inverse image of a 2n-th 
tensor in the Fock space denoted : (4>4>) n :. Note that : (00) := (00)— 2(7 These 
variables are orthogonal in I? . Set l x = l x — o x be the centered occupation 
field. Note that an equivalent formulation of proposition 5 is that the fields 
| : 00 : and L\ have the same law. 

Let us now consider the relation of higher Wick powers with self intersec- 
tion local times. 

9 Decompositions 

If D C X and we set F = D c , the orthogonal decomposition of the Dirichlet 
norm e(f) into e D (f — H F /) + e(H F f) (cf [6] and references) leads to the 
decomposition of the Gaussian field mentionned above and also to a decom- 
position of the Markov chain into the Markov chain killed at the exit of D 
and the trace of the Markov chain on F. 

Proposition 11 The trace of the Markov chain on F is defined by the 
Dirichlet norm e^(f) = e(H F f) , for which 

^x^y = Cx,y + ^2 Cx,aCb,y[G D ] a,b 
a,beD 

^ { x F} = ^- E c ^C b , x [G D ] a ' b 

a,beD 

and 

Z e = Z e DZ £ {F} 

Proof: The first assertion is well known. For the second, note first that 
for any y E F, [H F ]y = l x=y + l D (x) T,^ D [G D \ x ' b C b , y . Moreover, e(H F f) = 
(/, H F '/) e and therefore 

\ { / } = e {F} (l {x} ) = e(l {x} ,H F l {x} ) = X x - J2aeD C ^A HF }x = K ~ 

Then for distinct x and y in F, 

C x F y } = - (l {x} , l{y}) elF} = -(l {x} ,H F l {y} ) e = C x , y + J2a C xA HF \y = 
Cx,y+ Y.a,beDCx,aCb,y[G ]°' • 

Finally, note also that is the restriction of G to F. as for all x,y £ F, 
(G'' y , l{x}) e { F } — (G'' y , [H F ] x ) e = l{ x =y}- Hence the determinant decompo- 
sition already used in yields the final formula. 
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The cases where F has one point was already treated in section 3-2. 
The transition matrix [P^]^ can also be computed directly and equals 
Py+ EaMD P : P SV DU{x} ] a b = p y+ EaMD P :C b JG D ^r> b . The calcu- 

C {F} 

lation of-^j- yields a decomposition in two parts according whether the jump 
to y occurs from x or from D. 

If we set e x = e + IHI^^ and denote [e x ]^ by e^- F ' x ^ we have 

a,b 

and 

a,b 

More generally, if e* is such that C* = C on F x F, and A = A* on F 
we have: 

c*P = c x , y + J2cLMG' x D ] a ' b 

a,b 

and 

^ {F} = K-Y, c ^a, x [G' x D r) 

a,b 

A loop in X which hits F can be decomposed into a loop in F and 
its excursions in D which may come back to their starting point. 

Set v x,y = Cx,y5® + J2a,beD C x ,aCb,yfl>D aIld V D = K^'^+Y^l K ™ (J2a,beD C x , a C b>x ^ ) 

Here [i£ denotes the bridge measure (with mass [G D ] a,fe associated with e D . 
Note that ^^(1) = Ci F y and v^iX) = 'jpT- We get a decomposition of /i 

into its restriction yP to loops in D (associated to the process killed at the 
exit of D), a loop measure defined on loops of F by the trace of the 
Markov chain on F, measures on excursions in D indexed by pairs of 
points in F and measures on finite sequences of excursions in D indexed 
by points of F . 

Conversely, a loop of points in F (possibly reduced to a point), a 
family of excursions 7£ i ,£ i+1 attached to the jumps of and systems of i.i.d. 
excursions , 7^ attached to the points of defines a loop A(ft F \ (7&,£ i+1 ), (7&))- 
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Note excursions can be empty. Then fi — pP is the image measure by A of 
(J, {F} (dl {F} ) n(^ i+1 )( d 7&,&+i) Il^( d 7Ci), denoting v the normalised mea- 
sures jfa. 

The Poisson process = {l^ F \l G £ Q } has intensity and is 

independent of £®- 

In particular, if x is a measure carried by D, we have: 



More generally 

Proposition 12 if C* = C on F x F , and A = A* on F 

E£-(IL^FxF[^] JV -e-^ C ^ # - A -)|d J?> ) = 

The proof can be done by decomposing all e' into e* + (e' — e#), with 
(e' — e # ) carried byF 

These decomposition formulas extend to include a current c<j provided it 
is closed (i.e. vanish on every loop) in D. In particular, it allows to define 
uj f such that: 

Z e ,u> = Z e DZ e {F}^F 



10 Reflection positivity and Hilbert space 

Let us fix a. In view of physical applications, it is appropriate to assume 
that X is the union of two parts X ± exchanged by an involution p under 
which e is invariant. Each configuration C a of loops induces a configuration 
A of loops in X° = X + n X~. Given a function F on loops configuration in 
X + , it follows from the previous proposition the following 

Corollary 13 E(F(C a \ x+ )\A) = E(F o p(C a \ x -)\A) 
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so that the reflection positivity (also called physical positivity) property 
holds: 

||F|| = E(F(Fop)) > 0. 

The physical space is the quotient space modulo functionals of zero norm. 
It identifies with L 2 functionals of A. Osterwalder-Schrader-type construc- 
tion can be used to produce non commuting field observables. More precisely, 
after extending the framework to infinite spaces (see section below), one can 
assume for exemple X has a product structure S x Z and that the time 
translation r and the time reversal p leave e invariant. Then r induces a 
self adjoint contraction T of the physical space, hence a Hamiltonian log(T) 
and by complex exponentiation, a unitary dynamic U. Non commuting ob- 
servables are obtained by conjugation of an observable by the opeators U n . 
This extends the construction of the relativistic non commuting quantum 
free field observables out of the Euclidean Gaussian field. 

11 The case of general Markov processes 

We now explain briefly how some of the above results will be extended to 
a symmetric Markov process on an infinite space X. The construction of 
the loop measure as well as a lot of computations can be performed quite 
generally, using Dirichlet space theory Let us consider more closely the 
occupation field I. The extension is rather straightforward when points are 
not polar. We can start with a Dirichlet space of continuous functions and 
a measure m such that there is a mass gap. Let P t the associated Feller 
semigroup. Then the Green function is well defined as the mutual energy of 
the Dirac measures 5 X and 5 y which have finite energy. It is the covariance 
function of a Gaussian Markov field <p x , which will be associated to the field 
/ of local times. of the Poisson process of random loops whose intensity is 
given by the loop measure defined by the semigroup P t . More precisely, 
Propositions 1 and 5 still hold (x being defined as a Radon measure with 
compact support on X) as long as the continuous Green function G will be 
locally trace class. This will apply to exemples related to one dimensional 
Brownian motion or to Markov chains on countable spaces. 

When points are polar, one needs to be more careful. We will con- 
sider only the case of the two and three dimensional Brownian motion in 
a bounded domain killed at the boundary, i.e. associated with the classical 
energy with Dirichlet boundary condition. The Green function is not locally 
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trace class but it is still Hilbert-Schmidt which allows to define renormalized 
determinants det2 (Cf [12]) and to extend the statement of proposition 5 
to the centered occupation field and the Wick square : 00 : of the general- 
ized Gaussian Markov field 0. These three generalized fields are not defined 
pointwise but have to be smeared by measures of finite energy x such that 
J G x ' y x(dx)x(dy) < oo. The centered occupation field I is defined as follows: 
Let Af be the additive functional associated with x of finite energy. Then 

(j>iXj 1S defined as lim £ ^ Jq T ~ £ ^ dAf — (J>o(fo T ~^ dAf) which converges in 

L 2 (/z ). It is an intrinsic quantity. We then have 

Proposition 14 a) The centered occupation field C\ and the Wick square 
\ : 00 : have the same distribution. 

b) E(e-(^ x )) = det 2 (G x G- 1 ) a 

To justify the use of det 2 , note that in the finite case det(J + GM X ) = 
det 2 (J + GM x )e x ^ where we recall that a x 1 is the capacity of x, which 
vanishes now since x is polar. 

In two dimensions, higher Wick powers of 00 are associated with self 
intersection local times of the loops. 

Let us now consider currents. We will restrict our attention to the two 
dimensional Brownian case, X being an open subset of the plane. Currents 
can be defined by divergence free vector fields, with compact support. Then 
Jj^uj and J x (09 w — (pd^fydx are well defined square integrable variables 
(it can be checked easily in the case of the square by Fourier series). The 
distribution of the centered occupation field of the loop process "twisted" 
by the complex exponential exp(^ g£ Ii^ 00 + ^(IMI 2 )) appears to be the 
same as the distribution of : 00 : "twisted" by the complex exponential 
exp( J x (0^0- <t>dJ4>)dx) (Cf[10]). 

These points, among others, will be developped in a forthcoming article. 
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